A=list.files(path="C:/本科生教学/大三/R/data for class3/GSE67835",pattern = "csv")#读取文件夹中拓展名为.csv
setwd("C:/本科生教学/大三/R/data for class3/GSE67835")#设置路径
n = length(A) 
B = read.table(file = A[1]) #读取A中的第一个
for (i in 2:n)#从第二个开始循环
{
  C= read.table(file = A[i])
  B = merge(B,C,by.x=1,by.y=1)
}
list=substring(A,1,10)#取A文件名中的1到10长度的文件名编号
list=c("gene",list)#因为B中有一列为基因名字，所以给它加上列名，防止后面文件编号放在基因的列名上
colnames(B)=list#把B的列名进行更改
rownames(B)=B[,1]#把B的第一列作为行名
B=B[,-1]#删除第一列

E=read.table("C:/本科生教学/大三/R/data for class3/GSE67835/GSE67835-GPL15520_series_matrix.txt",skip = 36,nrows = 15)
#跳过36行
F=E[9:11,]
B=t(B)#把B的行名改为列名，列名改为行名
colnames(F)=list#把F的列名进行更改
F=F[,-1]#删除F第一列
G=rbind(B,F)
write.csv(G,file="实践4.csv")


H=read.csv("C:/本科生教学/大三/R/实践3/ADdate实践3-2.csv")
H=H[,-1]
rownames(H)=H[,1]
H=H[,-1]#必须将H改为这种形式
I1=apply(H, 1,function(x){sum(x==0)})
I2=apply(H, 2,function(x){sum(x==0)})#必须==，这是附属数值的意思
I1
I2

